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ABSTRACT 

Open clusters potentially provide an ideal environment for the search for transiting extrasolar planets 
since they feature a relatively large number of stars of the same known age and metallicity at the 
same distance. With this motivation, over a dozen open clusters are now being monitored by four 
different groups. We review the motivations and challenges for open cluster transit surveys for short- 
period giant planets. Our photometric monitoring survey (EXPLORE/OC) of Galactic southern open 
clusters was designed with the goals of maximizing the chance of finding and characterizing planets, 
and of providing for a statistically valuable astrophysical result in the case of no detections. We use 
the EXPLORE/OC data from two open clusters NGC 2660 and NGC 6208 to illustrate some of the 
largely unrecognized issues facing open cluster surveys including severe contamination by Galactic field 
stars (> 80%) and relatively low number of cluster members for which high precision photometry can 
be obtained. We discuss how a careful selection of open cluster targets under a wide range of criteria 
such as cluster richness, observability, distance, and age can meet the challenges, maximizing chances to 
detect planet transits. In addition, we present the EXPLORE/OC observing strategy to optimize planet 
detection which includes high-cadence observing and continuously observing individual clusters rather 
than alternating between targets. 

Subject headings: techniques: photometric, surveys, planetary systems, open clusters and associations: 
general, individual (NGC 2660, NGC 6208) 



1. INTRODUCTION 

The EXPLORE Project (Mallen-Ornelas et al. 2003; Yee 
et al. 2003) is one of currently about 20 ongoing surveys 1 
with the aim to detect transiting close-in extrasolar giant 
planets (CEGPs; also referred to as '51-Peg type' or 'hot 
Jupiters', i.e., planets with radius of order Jupiter-radius, 
orbital period of one to a few days, and transit durations of 
a few hours) around Galactic main-sequence stars. Transit 
studies explore a different parameter space in the search for 
extrasolar planets from the very successful radial- velocity 
or "wobble" method. This method based on detecting 
planets via the radial motions of their parent star caused 



by the star's motion about the common center of mass 
(see for example table 3 in Butler et al. 2002). Fainter 
(and thus more) stars can be monitored photometrically 
than spectroscopically. Thus, more distant environments 
can be probed for the existence of extrasolar planets. 

All transiting planets have a measured radius, based on 
transit depth and stellar radius. Knowledge of the planet's 
radius and mass plays an important role in modeling in- 
ternal structure of planets and hence the formation, evolu- 
tion, and migration of planetary systems (see for example 
Burrows et al. 2000; Guillot & Showman 2002; Baraffe 
et al. 2003, and references therein). 

Transiting planets are currently the only planets whose 



1 See http://star-www.st-and.ac.uk/~kdhl/transits/table.html, maintained by K. Home. 
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physical characteristics can be measured. In addition to 
mass and radius, several parameters can be constrained 
from follow-up measurements. For example, the fact that 
a transiting planet will be superimposed on its parent star 
can be used to determine constituents of the planet's at- 
mosphere by means of transmission spectroscopy (Char- 
bonneau et al. 2002; Vidal-Madjar et al. 2003), or put 
constraints on the existence of planetary moons or rings 
(e.g., Brown et al. 2001). In addition, the secondary eclipse 
can provide information about the planetary temperature 
or its emission spectrum (Richardson et al. 2003a, b). 

At the time of writing, six transiting planets are known. 
HD209458b (Charbonneau et al. 2000; Henry et al. 2000; 
Brown et al. 2001) was discovered by radial-velocity mea- 
surements (Henry et al. 2000; Mazeh et al. 2000) and the 
transits were discovered by photometric follow-up. OGLE- 
TR-56 (Udalski et al. 2002a,b; Konacki et al. 2003a) was 
the first planet discovered by the transit method, and the 
first of currently four planets based on photometry of the 
OGLE-III survey (Udalski et al. 2002a,b, 2003): OGLE- 
TR-113 (Bouchy et al. 2004; Konacki et al. 2004), OGLE- 
TR-132 (Bouchy et al. 2004), and OGLE-TR-111 (Pont 
et al. 2004). Very recently, Alonso et al. (2004) found 
a further transiting planet, TrES-1, using telescopes with 
10cm apertures. Over 20 transit searches are currently 
ongoing to find more planets. 

As part of the EXPLORE 2 Project, we have recently be- 
gun a survey - EXPLORE/OC - of southern open clusters 
(OCs) with the aim of detecting planetary transits around 
cluster member stars. During the course of ^3 years, we 
hope to conduct searches of up to 10 OCs using the Las 
Campanas Observatory (LCO) lm Swope Telescope. To 
date we have monitored five OCs (see §6). 

In addition to EXPLORE/OC, at the time of writing 
there are currently three OC planet-transit surveys under- 
way 3 : 



The OC NGC 6791 was furthermore monitored for plan- 
etary transits by Bruntt et al. (2003). As all these sur- 
veyed OCs are located the northern hemisphere, EX- 
PLORE/OC 7 is currently the only OC survey operating 
in the south where most of the Galactic OCs are located. 

With the growing number of open star cluster sur- 
veys this publication describes incentives, difficulties, and 
strategies for open cluster planet transit surveys, thereby 
including a discussion on transit surveys in general. We 
use data from the first two targets (NGC 2660 and NGC 
6208) from our program to illustrate the major issues for 
OC transit surveys. 

The concept and advantages of monitoring OCs for the 
existence of transiting planets were originally described in 
Janes (1996). Written before the hot Jupiter planets were 
discovered, Janes (1996) focused on 12-year period orbits 
and long-term photometric precision required to determine 
or put useful limits on the Jupiter-like planet frequency. 
This paper is intended to be a modern version of Janes 
(1996) based on the existence of short-period planets and 
practical experience we have gained from both the EX- 
PLORE and the EXPLORE/OC planet transit surveys. 
Section 2 motivates OC transit surveys. Section 3 ad- 
dresses the challenges facing transit surveys in general, 
and §4 addresses challenges specifically facing OC tran- 
sit surveys. Section 5 focuses on strategies to select OCs 
which are most suited for transit surveys and which min- 
imize challenges described in the previous sections. The 
EXPLORE/OC strategies concerning target selection, ob- 
serving methods, photometric data reduction, and spec- 
troscopic follow-up observations are described in §6, §7, 
and §8, respectively. These Sections contain relevant pre- 
liminary results on EXPLORE/OC's first two observed 
clusters NGC 2660 and NGC 6208. We summarize and 
conclude in §9. 



• The Planets In Stellar Clusters Extensive Search 
(PISCES 4 ) reported the discovery of 47 and 57 low- 
amplitude variables in the open clusters NGC 6791 
(Mochejska et al. 2002) and NGC 2158 (Mochejska 
et al. 2004), respectively. 

• The University of St. Andrews Planet Search 
(UStAPS 5 ) has monitored the OCs NGC 6819 
(Street et al. 2002) and NGC 7789 (Bramich et al. 
2003) and published data on variable stars in NGC 
6819 (Street et al. 2003). 

• The Survey For Transiting Extrasolar Planets In 
Stellar Systems (STEPSS 6 ), described in Burke 
et al. (2003) and Gaudi et al. (2002). They have an- 
alyzed monitoring data of the OC NGC 1245 (Burke 
et al., in preparation) and determined its fundamen- 
tal parameters (Burke et al. 2004) . Analysis of their 
data on NGC 2099 and M67 is currently ongoing. 

2 http://www.ciw.edu/seager/EXPLORE/explore.htm 

3 see also http://www.ciw.edu/kaspar/OC_transits/OC_transits.html 

4 http:/ /cfa-www. harvard.edu/~bmochcjska/PISCES 

5 http://crux.st-and.ac.uk/~kdhl/ustaps.html and http://star-www.st- 

6 http://www.astronomy.ohio-state.edu/~cjburkc/STEPSS 

7 http://www.ciw.edu/seager/EXPLORE/opcn_clusters_survey.html 



2. MOTIVATION FOR OPEN CLUSTER PLANET TRANSIT 
SEARCHES 

Open clusters present themselves as "laboratories" 
within which the effects of age, environment, and espe- 
cially metallicity on planet frequency may be examined. 
Evidence that planet formation and migration are corre- 
lated with metallicity comes from radial velocity planet 
searches (Fischer & Valenti 2003). The fact that no plan- 
etary transits were discovered in the monitoring study of 
47 Tuc by Gilliland et al. (2000) may be due to its low 
metallicity, or alternatively due to the high-density envi- 
ronment in a system such as a globular cluster (or due to 
both effects). The less-crowded OCs of the Milky Way of- 
fer a range of metallicities and may thus be further used 
to disentangle the effects of metallicity versus high-density 
environment upon planet frequency. 

Monitoring OCs for the existence of planetary transits 
offers the following incentives (see also Janes 1996; Char- 
bonneau 2003; Lee et al. 2004; von Braun et al. 2004): 
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1. Metallicity, age, distance, and foreground reddening 
are either known or may be determined for clus- 
ter members (more easily and accurately than for 
random field stars; see §5 and, e.g., Burke et al. 
2004). Thus, planets detected around open cluster 
stars will immediately represent data points for any 
statistic correlating planet frequency with age, stel- 
lar environment, or metallicity of the parent star. 

2. The planet-formation and planet-migration pro- 
cesses, and hence planet frequencies, may differ be- 
tween the OC, globular cluster, and Galactic field 
populations. Planet transit searches in OCs, to- 
gether with many ongoing transit field searches and 
GC surveys (e.g., the ground-based work on 47 Tuc 
by Weldrake et al. 2003, 2004), enable comparison 
between these different environments. 

3. Specific masses and radii for cluster stars may be 
targeted (within certain limits of other survey de- 
sign choices, see §7.2) in the planet-search by the 
choice of cluster distance and by adjusting expo- 
sure times for the target. 

3. MAIN CHALLENGES FOR TRANSIT SURVEYS 

Open cluster planet transit surveys are a subset of planet 
transit surveys and therefore have some important chal- 
lenges in common. Articulating these challenges is crucial 
in light of the fact that over 20 planet transit surveys have 
been operating for a few years (Home 2003), with only six 
known transiting planets of which five were discovered by 
transits. 

The most basic goals of any transit survey are to (1) 
detect planets and provide their characteristics, and (2) 
to provide (even in the case of zero detections) statistics 
concerning planet frequencies as a function of the astro- 
physical properties of the surveyed environment. Although 
the most important considerations for designing a success- 
ful transit survey were presented in Mallen-Ornclas et al. 
(2003, M03 hereafter) for the EXPLORE Project, we sum- 
marize and provide updates to the three key issues: num- 
ber of stars, detection probability, and blending. For OC 
surveys, these issues (with the exception of blending) can 
be optimized or overcome by a careful survey strategy, par- 
ticularly by the selection of the target OC (see §5 through 
7). 

3.f . Maximizing Number of Stars with High Photometric 

Precision 

Any survey's goal should be to maximize the number of 
stars for which it is possible to detect a transiting planet. 
We discuss the three most important aspects below: the 
astrophysical frequency of detectable transiting planets, 
the probability with which existing planetary transits are 
observed, and the number of stars for which the relative 
photometric precision is sufficiently high. 

The frequency of detectable transiting planets is calcu- 
lated by considering the astrophysical factors: frequency 
of CEGPs around the surveyed stars; likelihood of the ge- 
ometrical alignment between star and planet necessary to 
detect transits; and binary fraction. We assume a planet 
frequency around isolated stars of 0.7% for planets with 



semi-major axis a ~ 0.05 (Marcy et al. 2004; Naef et al. 
2004). Of those CEGP systems, approximately 10-20% 
(probability <~ R^/a) would, by chance, have a favorable 
orientation such that a transit is visible from Earth. We 
assume that planets can only be detected around single 
stars and conservatively adopt a binary fraction of 50%; 
although there are known planets orbiting binary stars and 
multiple star systems (Zucker & Mazeh 2002; Eggenberger 
et al. 2004), their detection by transits would be difficult 
due to a reduced photometric signature in the presence 
of the additional star. Combining the above estimates, 
we arrive at the at the value of 1 star in 3000 (a ~ 0.05 
AU) having a hot Jupiter planetary transit around a main 
sequence star. 

The probability of detecting an existing transiting hot 
Jupiter (1/3000) applies only to stars for which it is possi- 
ble to detect a planetary transit given the observational 
setup of the survey. This number of "suitable" stars 
is frequently equated to the number of stars with high 
enough relative photometric precision to detect the tran- 
siting planet (see Fig. 1). Planet transit surveys in gen- 
eral reach photometric precision sufficiently high to de- 
tect Jupiter-sized transiting planets around main sequence 
stars (see Fig. 1) for up to 40% of stars in their survey, de- 
pending on crowdedness and other factors. For example, 
the EXPLORE search reached relative photometric preci- 
sion better than 1% on 37,000 stars from 14.5 < I < 18.2 
out of 350,000 stars down to / = 21 (M03); OGLE-III 
reached better than 1.5% relative photometric precision 
on 52,000 stars out of a total of 5 million monitored stars 
(Udalski et al. 2002b); the Sleuth survey (O 'Donovan et al. 
2004), reaches better than 1.5% relative photometry over 
the entire dataset on the brightest 4000 stars out of 10,000 
in their 6x6 degree 2 field, using an automated 10 cm tele- 
scope with a 6 degree square field of view. 

The real number of stars suitable for planet transit de- 
tection, however, is not equivalent to the number of stars 
with 1% relative photometry. One may see from Fig. 1, 
for instance, that a Jupiter-sized planet would cause a 2% 
eclipse around a parent MO-star. Furthermore, Pepper & 
Gaudi (in preparation) find that, if a planet with given 
properties around a cluster member star on the main se- 
quence produces a detectable transit signature, a planet 
with identical properties orbiting any other main-sequence 
cluster member will produce a detection of approximately 
the same signal-to-noise-ratio (SNR) unless the sky flux 
within a seeing disk exceeds the flux of the star. Since 
most transit surveys aim to find planets of approximately 
Jupiter-size around stars whose radii are close to, or less 
than, a solar radius, the number of stars with 1% relative 
photometric precision can therefore be regarded as a lower 
limit to the number of stars suitable for transit detection. 
For the rest of this publication, we thus use this number 
as a proxy for the number of stars around which we (or 
other transit surveys) can detect planets. 

3.2. Probability to Detect an Existing Transit 

The actual observed hot-Jupiter transit frequency will 
be lower than 1/3000 due to the probability with which 
an existing transit would be observed two or more times 
during an observing run. This probability, which we call 
P V i S , is equivalent to the window function of the obser- 
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vations. Although the probability function has been de- 
scribed in detail before (Borucki & Summers 1984; Gaudi 
2000, and M03), we extend the discussion to include the 
recently discovered class of 1-day period planets, as well 
as considering different metrics for probability to detect 
an existing planet. In all of our P vis simulations we as- 
sume the simplified case of a solar mass, solar radius star 
with the planet crossing the star center, thus focusing on 
stars we are most interested in. The transit duration is 
then related to the planet's period by ^duration = PRq/tto, 
(typically a few hours for a period of a few days). 

In panel a of Figure 2 we show the P vis for detecting 
existing transiting planets with different orbital periods 
under the requirement that two or more full transits must 
be observed. We consider different runs (7, 14, 21 nights) 
of consecutive nights with 10.8 hours of uninterrupted 
observing each night with 5-min time-sampling. This is 
equivalent to approximately 125 observations of a given 
cluster per night. The P V i S for 1- to 2-day period planets 
is basically complete for the 14- and 21-night runs, while 
the P v is is markedly lower for a 14-night run compared to 
a 21-night run for planets with periods between 2 and 4 
days. 

We now turn to P vis for a transit detection strategy 
where it is not necessary to detect a transit in its entirety 
during a single observing night. Instead, the strategy re- 
quires a transit to be detected in phased data (from at 
least two individual transit events). Such a P V i S (which we 
call PviSphased) is relevant for transit detections based on 
period-folding transit-searching algorithms (for instance 
with data covering partial nights or a strategy of alter- 
nating targets throughout the course of a night). With 
a period-folding algorithm, each individual transit need 
not be fully sampled. In order to quantify Pvis p hased, we 
specify that the phased transit must be sampled by at least 
N points. Since a typical duty cycle of a transit is on the 
order of a few percent, we choose N = 20, 40, or 60 to rep- 
resent light curves with a total of a few hundred to a thou- 
sand data points (the phased OGLE planets' light curves 
typically have a few tens of data points obtained during 
transit). Pvis p hased is then calculated to be likelihood (as 
a function of period) that at least N in-transit points are 
accumulated for observing runs of different lengths and 
different observing cadences. 

Note that in reality a detection of a planet transit de- 
pends on the number of photons observed during the tran- 
siting phase. This number of photons is contained in the 
combination of the SNR per individual data point and the 
number of data points (during any transit). A back-of-the- 
envelope calculation would give a transit SNR for a Am <~ 
2%-dcpth transit with M = 20 data points and a relative 
photometry precision of rms <~1%: 

5M,^x^ = V20x^~9. (1) 
rms 0.01 

For comparison with the two- full-transit P V i S , we show 
in Fig. 3 Pvis p h aS ed for N = 20,40,60 by a solid, dotted, 
and dashed line, respectively The four panels represent 
different observing strategies: 

• Panel a shows Pvis p hased for 21 nights (10. 8h) with 
5-minute time-sampling, resulting in a light curve 
with around 2700 data points. If the data SNR is 



high enough for 20 data points during transit to 
constitute a detection, then Pvis p h aS ed is high for 
all transit periods between one and five days. If, 
in contrast, the data SNR is lower, and 40 or 60 
points per transit are required for detection, then 
Pvisphased is low for P > 2 days. 

• Alternatively, one could imagine a strategy of alter- 
nating between cluster fields (to increase the num- 
ber of monitored stars), in which case the observing 
cadence is reduced. Panel b shows Pvis p hased when 
observing for 21 nights with a 1 5-minute cadence (~ 
900 measurements in the light curve). The proba- 
bility to detect transits with N = 20 is very low for 
P > 2 days, and it is zero for TV = 40 or 60. 

• Panel c shows Pvis phase d for 40 nights of contin- 
uous observing (10. 8h) with an observing cadence 
of 5 minutes (~ 5200 data points). For N = 20 
and N = 40, Pvis p hased is close to complete for all 
shown periods. 

• Panel d shows Pvis p hased for 40 nights of observing 
with a 15-minute cadence (again simulating a strat- 
egy of alternating between cluster fields; ~ 1700 
data points). Pvis p hased is very low for N > 20, in- 
dicating that, in order to be able to observe with a 
15-cadence, many more than 40 nights are needed if 
more than 20 data points are required for a transit 
detection. 

We conclude that the ability to particularly detect 
longer-period (P > 2 days) planets depends on observ- 
ing strategy. For the rest of this paper we adopt the P vis 
criterion of seeing two full transits which is a good strategy 
for a limited number (~ 20) of observing nights. 

3.3. Blending and False Positives 

Blending in a planet transit light curve due to the pres- 
ence of an additional star is a serious challenge inherent 
in planet transit surveys, one that has only recently been 
gaining recognition. If the light of multiple stars arc in- 
terpreted as being due to one individual star, then the 
relative depth of any eclipse will be decreased. This "light 
pollution" may either cause (1) an eclipsing binary sys- 
tems to mimic a more shallow transiting planet signal, or 
(2) a true planet's transit signal's depth to be decreased to 
a fraction of its already very small amplitude, rendering it 
harder or even impossible to detect. 

Blending can be caused either by optical projection in 
crowded fields, or by physically associated stellar systems. 
The crowdedness can be considered in the choice of target 
and observational setup. Blending in spatially unresolved, 
physically associated systems generally consists of a wide 
binary of which one component hosts an additional close- 
by stellar companion or a transiting planet. The compo- 
nent without the close-by stellar or planetary companion 
would produce the polluting light. Although we have ac- 
counted for binary stars in our probability estimate (§3.1), 
the contribution due to this kind of "false positive" may 
be larger due to the unknown wide-binary component dis- 
tances. 
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Recently the effect of blending on the probability of de- 
tecting planets has been addressed by several authors (e.g., 
M03; Brown 2003; O'Donovan ct al. 2004; Konacki et al. 
2003b) in the context of causing false positives. Several so- 
lutions have been proposed to avoid false positive transit 
candidates that are actually blended star light curves. Sea- 
ger & Mallcn-Ornclas (2003) show that one may eliminate 
some false positives due to blending with photometric data 
alone if the light curve is of sufficient relative photomet- 
ric precision and the observing cadence is high enough to 
clearly resolve the individual temporal components of the 
transit. Using spectroscopic data, other solutions include 
a careful modeling of the additional star properties to de- 
tect a second cross correlation peak caused by a physically 
associated star (M03; Konacki et al. 2003b; Kotredes et al. 
2003; Torres et al. 2004b). Finally, estimates for blending 
(associated either with chance alignment of foreground or 
background stars, or physical triplets) effects on the prob- 
ability of detecting existing transits can be quantified for 
individual surveys, as done by Brown (2003) for shallow 
wide-field transit surveys. 

4. MAIN CHALLENGES FOR OPEN CLUSTER TRANSIT 
SURVEYS 

The difficulties and challenges involved in searching for 
planetary transits specifically in OCs are: the fixed and 
somewhat low number of stars in an open cluster, deter- 
mining OC cluster membership in the presence of signif- 
icant contamination, and differential reddening along the 
cluster field and along the line of sight. We outline these 
aspects individually below. 

1. The Number of Monitored Stars The num- 
ber of monitored stars is typically lower than in 
rich Galactic fields (in part due to the smaller field 
size of the used detectors), reducing the statisti- 
cal chance of detecting planets. Open clusters can 
have up to <~ 10,000 member stars (Friel 1995), de- 
pending on the magnitude range taken into con- 
sideration. Only a subset of these stars, perhaps 
10-20%, however, will be observed with sufficient 
relative photometric precision to detect transits (cf. 
§3.1). The number of these stars in rich OCs is com- 
parable to the number in wide-field, shallow transit 
surveys, e.g., Sleuth: ~4000 stars in 6x6 square 
degrees of 9 < R < 16 (O'Donovan ct al. 2004); 
WASPO: < 3000 stars in 9x9 square degrees of 
broadband magnitude between 8 to 14 (Kane et al. 
2004). The richest deep galactic fields surveyed 
have many more stars with high relative photomet- 
ric precision (§3.1). 

2. Cluster Contamination Determining cluster 
membership of stars in the OC fields without spec- 
troscopic data or proper motion information is dif- 
ficult due to significant contamination by Galac- 
tic field stars, since the clusters are typically con- 
centrated toward the Galactic disk. For example, 
Street ct al. (2003) estimate the contamination of 
Galactic field stars in their study of NGC 6819 to 
be around 94%. A study by Nilakshi et al. (2002) 
calculated the average contamination in the fields of 
38 rich OCs to be 35% in the inner regions and 80% 



in the "coronae" of the clusters. Furthermore, if the 
target is located such that the line of sight includes 
a long path through the Galaxy (e.g., low Galactic 
latitude and longitude towards the Galactic bulge), 
background giants may start polluting the sample 
of stars with apparent magnitudes monitored with 
high relative photometric precision (see for exam- 
ple the discussion in Street et al. 2003). Getting 
a handle on the issue of contamination is vital for 
OC surveys since any statistical statements about 
the result will need to be based on estimates of sur- 
veyed cluster members. 

3. Differential Reddening Differential reddening 
across the cluster field and along the line of sight 
can make isochrone fitting (and subsequent deter- 
mination of age, distance, and mctallicity) difficult. 
Ranges of AEb-v ~ 0.2 or higher across fields-of- 
view of 10-20 arcmin on the side are not uncommon 
(see for example the studies by Munari & Carraro 
1996; Raboud et al. 1997; Rosvick & Balam 2002; 
Carraro & Munari 2004; Villanova et al. 2004; Pris- 
inzano et al. 2004). For an Ry = 3.1 reddening-law, 
a differential AEb-v ~ 0.2 corresponds to a differ- 
ential AV ~ 0.6 and Al ~ 0.3 (Cardelli et al. 1989). 
The calculated effective temperature of a solar- 
metallicity main-sequence star with V — I ~ 0.8 
would vary by about 500 K for a differential red- 
dening effect of AEb-v ~ 0.2 (Houdashelt et al. 
2000). It should, however, also be noted that some 
OCs do not seem to suffer from differential redden- 
ing, such as NGC 1245 as examined in Burke et al. 
(2004) and NGC 2660 in our preliminary analysis 
of its CMD. 

5. OPEN CLUSTER SELECTION 

Open cluster target selection can help overcome or re- 
duce some of the main challenges of OC planet transit 
surveys described in §3 and §4. More specifically, careful 
cluster selection can help maximize the number of stars, 
maximize the probability of detecting existing transits, 
and reduce line-of-sight and differential reddening. Most 
importantly, cluster selection allows for targeting specific 
spectral type for a given telescope and observing cadence. 

The biggest challenge in the selection of target clusters 
is the paucity of data about many OCs. The physical 
parameters of the cluster, such as distance, foreground 
reddening, age, and metallicity are frequently either not 
determined, or there exist large uncertainties in the pub- 
lished values. For example, out of approximately 1100 
associations of stars designated as OCs, many only have 
identified coordinates, approximately half have an estab- 
lished distance, and about 30% have an assigned metal- 
licity (WEBDA database Mermilliod 1996). Additional 
difficulties arise when independent studies arrive at differ- 
ent values for any of the parameters. For instance, the 
metallicity of NGC 2660 was determined to be as low as 
-1.05, and as high as +0.2 (see discussion in the introduc- 
tion of Sandrelli ct al. 1999). The very important criterion 
of richness tends to be even less explored than the other 
physical parameters, probably due to the significant con- 
tamination from field stars that these OCs tend to suffer. 
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of sight (Schlegel et al. 1998). Differential reddening will 
cause the main sequence of the cluster to appear broad- 
ened (e.g., von Braun & Mateo 2001). This, in turn, will 
cause large errors in the determination of cluster parame- 
ters such as age, metallicity, etc, by means of isochronc fit- 
ting. Furthermore, a broad main sequence will make any 
attempts to estimate contamination (based on isochrone 
fitting; see for example Mighell et al. 1998; von Hippel 
et al. 2002) more challenging. 

5.3. Cluster Age 

The consideration of cluster age in the OC selection 
is not as crucial as richness, observability, and distance. 
However, choosing both younger and older OCs may im- 
pose different challenges with respect to the transit-finding 
process. 

Stellar surface activity, which would introduce noise into 
the light curve of a given star, decreases with age. As they 
age, stars lose angular momentum and thus magnetic ac- 
tivity on their surfaces (see for example Donahue 1998; 
Wright 2004, and references therein). The photometric 
variability for a sample of Hyades OC stars was found to 
be on the order 0.5-1% (Paulson et al. 2004) with periods 
in the 8-10 day range (the Hyades cluster has an age of 
~650 Myr; see Perryman et al. 1998; Lebreton et al. 2001). 
While these photometric variations do not necessarily rep- 
resent a source of contamination in the sense of creating 
false positives, they nevertheless will introduce noise into 
the stellar light curves and thus render existing transits 
more difficult to detect. 

Stellar surface activity is potentially an issue not just 
for host star photometric variability, but more so for ra- 
dial velocity follow-up. Paulson et al. (2004) found that 
radial velocity rms due to rotational modulation of stellar 
surface features can be as high as 50 m/s and is on aver- 
age 16 m/s for the same sample of Hyades stars. Further 
such correlation between radial velocity rms and photo- 
metric variability was found by Queloz et al. (2001) who 
observed a v r amplitude of ~ 180 m/s for HD 166435 (age 
200 Myr), a star without a planet. The associated pho- 
tometric variability with a period of around 3.8 days is of 
order 5%. 

This variability issue favors older OCs as targets, partic- 
ularly since most of the decrease in surface activity occurs 
between stellar ages between 0.6 Gyrs and 1.5 Gyrs (Pace 
& Pasquini 2004). A 16 m/s rms may not be a problem 
for deep OC surveys; for short-period Jupiter-mass planets 
relatively large radial velocity signatures are expected and 
the faint stellar magnitudes limit radial velocity precision 
to ~ 50-100 m/s (Konacki et al. 2003a, 2004; Bouchy et al. 
2004) using currently available telescopes and instrumen- 
tation. 

Older star clusters offer an additional advantage: in 
general older OCs are richer and more concentrated, and 
therefore offer a larger number of member stars to be sur- 
veyed (Friel 1995). On the other hand, some old OCs ap- 
pear to be dynamically relaxed and mass-segregated (such 
as NGC 1245; Burke et al. 2004), and in the case of NGC 

8 http://obswww.unigc.ch/wcbda/ 

9 To optimize this observability for our potential cluster targets, one may use SKYCALC, written by J. Thorstensen and available at 
ftp://iraf.noao.edU/iraf/contrib/skycal.tar.Z 

10 Available, for instance, at http://cadcwww.hia.nrc.ca/cadcbin/gctdss 



In spite of the lack of OC data, a very useful place to 
start is the WEBDA 8 database (Mermilliod 1996). From 
the long list of potential OC monitoring targets, one can 
then start eliminating cluster candidates by applying the 
criteria we describe below. 

5.1. Cluster Richness and Observability 

Apart of its observability for a given observing run 9 , the 
most important selection criterion for a cluster is its rich- 
ness, simply to increase the statistical chance of detecting 
planets. The richness of the cluster field can be estimated 
by looking at sky-survey plots 10 of the appropriate region. 
Estimating the richness of the cluster itself is a much more 
difficult process since field star contamination is usually 
significant due to the typically low Galactic latitude of the 
Pop. I OCs (see below and Bramich et al. 2003; Street 
et al. 2003; Burke et al. 2004). One may use published 
cluster richness classifications, such as in Cox (2000), taken 
from Janes & Adlcr (1982). The WEBDA database also 
contains information for clusters in the 1987 Lynga catalog 
(online data published in Lynga 1995). The data published 
there, however, only present lower estimates for richness 
classes. The depths of the studies from which the rich- 
ness classes were derived may differ significantly from one 
study to the next. Thus, these classes are rough estimates 
only, and the best way to judge the richness of a cluster 
field is to rely on one's own test data obtained with the 
same setup as the one used for the monitoring study (see 
§7.5). 

5.2. Cluster Distance 

The distance to the target cluster is an important crite- 
rion for cluster selection for four reasons: (1) to ensure the 
cluster is sufficiently distant to fit into the field of view, (2) 
to allow RV follow-up of potential candidates, (3) to tar- 
get the desired range of spectral types for given observing 
conditions, and (4) to minimize reddening. 

Since all stars in an OC are at approximately the same 
distance, one may, with appropriate adjustment of the ex- 
posure time for given telescope parameters, cluster dis- 
tance, and foreground reddening, target certain spectral 
types of stars for high-precision photometry. We arc in- 
terested in GO or later spectral type/smaller stars, since 
early-type stars have larger radii which would make tran- 
sit detections more challenging. The transit depth in the 
light curve is, for small Am, simply given by 

AF ' ' • ■ - 



Am - 
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where Am and AF are the changes in magnitude and flux, 
respectively, and Fq is the out-of-transit flux of the parent 
star (see Fig. 1). In addition to featuring larger radii, early 
spectral types are fast rotators which exhibit broad spec- 
tral lines, making mass-determination of planetary com- 
panions more difficult. 

Foreground reddening, increasing with cluster distance, 
usually represents a proxy for the amount of differential 
reddening across the field of view as well as along the line 
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3680, for instance, evidence seems to point toward some 
resulting evaporation of low-mass stars over time (Nord- 
stroem et al. 1997). Since low-mass stars are the primary 
monitoring targets, some dynamically evolved OCs may 
actually be less favorable for observing campaigns. 

5.4. Other Criteria 

Given a sufficiently high remaining number of suitable 
open clusters after consideration of the previous four selec- 
tion criteria, cluster metallicity and galactic location are 
additional relevant selection criteria. 

• Range of metallicities. In order to be able to make 
quantitative statements about planet frequency as 
a function of metallicity of the parent star (§2), 
one needs to have a sample of clusters with varying 
metallicities. Surveys based on the radial-velocity 
method indicate that solar-neighborhood stars with 
higher metallicities are more likely to harbor plan- 
ets than metal-poor ones (Fischer & Valenti 2003). 
It may therefore be advantageous to favor highcr- 
metallicity clusters for monitoring studies to find 
planets. 

• The target's Galactic coordinates. On average, the 
closer the OC is to the Galactic disk, the higher the 
contamination due to Galactic field stars (see §3). 
Moreover, if the target is located close to, or even in 
front of, the Galactic bulge, contamination may be 
severe. It should be pointed out that background 
giants and subgiants will truly pollute the stellar 
sample since their radii are too large to detect plan- 
ets. Transit detections around main-sequence field 
stars at distances of less than, or roughly equal to, 
the OC distance are still possible and would be as 
scientifically valuable as a detection of a planetary 
transit as part of a dedicated field survey. 

6. explore/oc target-selection strategy 

EXPLORE/OC is a transit survey in open clusters op- 
erating with the LCO lm Swope Telescope with a field of 
view of 24 arcmin x 15 arcmin and a scale of 0.435 arc- 
sec/pixel. We have observed 5 open clusters to date: NGC 
2660 observed for ~ 15 nights in Feb 2003 (von Braun et al. 
2004); NGC 6208 observed for ~ 21 nights in May- June 
2003 (Lee et al. 2004); IC 2714 observed for - 21 nights 
in March/April 2004; NGC 5316 observed for - 19 nights 
in April 2004; and NGC 6253 observed for - 18 nights in 
June 2004 11 . 

Our /-band, high-cadence (~ 7 min, including 2 min 
readout time) photometric monitoring enables us to typ- 
ically attain 1% precision in our relative photometry for 
around 3000-5000 stars per cluster target field in the range 
14.5 < I < 17 (see Fig. 5). This number corresponds to 
a lower limit to the number of stars around which we can 
detect planetary transits (see §3.1). In the context of out- 
lining our survey strategies, we present some of our prelim- 
inary results of the studies of the open clusters NGC 2660 
and NGC 6208. In this Section we explain our approach 



to target selection, specifically designed to maximize the 
number of target stars of appropriate spectral type. 

6.1. Overall Potential Targets 

Our potential OC targets listed in Table 1 were cho- 
sen with the basic goal that we observe as many cluster 
member stars as possible at a sufficiently high photomet- 
ric precision and high-cadence of observations to detect 
CEGPs around them. Richness classes are given when- 
ever they were available. We note that we used these pub- 
lished richness classes only as a guideline (i.e., we gave 
extra considerations to OCs classified as rich, but did not 
necessarily discard any OCs classified as poor) and relied 
more on visual inspection and photometric analysis of sky 
survey images of the cluster regions. 

Targets in Table 1 were further selected based on the 
published estimates for distance and foreground redden- 
ing 12 . To select a cluster with a suitable distance we con- 
sider the preferred range of spectral types (G to M), our 
adopted relatively short exposure times (see §7.2), and the 
size of the LCO Swope telescope. 

As an example of how distance, exposure time, target 
spectral type, and reddening are related we use our OC 
NGC 6208. In Fig. 5, we show our photometric precision 
as a function of I magnitude of our NGC 6208 data (night 
15), obtained during May and June 2003 at the LCO lm 
Swope Telescope. We conservatively estimate that, with 
our exposure time of 300s per frame, we attain 1% preci- 
sion for a range of about 2.5 magnitudes (14.5 < I < 17). 
From the WEBDA database (see also Table 1) we find 
that the distance to NGC 6208 is 939 pc, and the fore- 
ground reddening is Eb-v = 0.210. Using the relation 
Ai = 1.94E B -v from Schlegel et al. (1998), we find that 
I = 17 for NGC 6208 cluster members corresponds to 
Mj = 6.73 which we call Mj limit in Table 1. Using table 
15.7 in Cox (2000), this corresponds to an MK spectral 
type of M0 or Ml. The bright limit, above which satura- 
tion will start to set in, would be at Mj ~ 4.2 which would 
correspond to an MK spectral type of approximately G5 
(Cox 2000). Once the range of spectral types for the mon- 
itored cluster members is determined, we can estimate the 
range of planetary radii which would be detectable (see 
Fig. 1). 



6.2. Potential Targets for a Given Observing Run 

For a given slot of observing time, we use Table 1 as 
the source from which we pre-select two or three potential 
observing targets for the run. The final target selection is 
then performed based on our own data, taken either dur- 
ing a previous observing run or at the very beginning of 
the observing run itself (see below). The main criterion at 
the pre-selection stage is the observability of the potential 
targets to maximize the time during which we can observe 
the OC. 

With all of the other constraints (distance, reddening, 
richness) on cluster selection, finding a cluster that is ob- 
servable all night long becomes challenging when observing 
runs are long. The main criterion for a successful transit 

11 We define here the number of nights as the number of at least partially useful nights during the monitoring campaigns. 

12 Note that our first target, NGC 2660, selected on the basis of its estimated richness and observability alone, turned out to have a relatively 
large distance and high foreground reddening. 
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search is maximizing the time during which we can ob- 
serve the respective target OC (to increase P V i S ), making 
clusters of numerically high southern declination prefer- 
able targets. 

6.3. Final Target Selection 

Our final target selection is based on the evaluation of 
VI test data (see Fig. 6) for the group of pre-selected 
clusters, involving the following steps: 

1. We create (at least roughly) calibrated VI CMDs 
of the potential target clusters based on our own 
test data. These data were obtained either during 
the beginning of the same observing run or during 
prior runs during photometric conditions and rea- 
sonable seeing, and have the same exposure time as 
for the eventual monitoring. Figure 6 shows these 
CMDs for the OCs NGC 6253 (left panel) and NGC 
6134 (right panel), both of which were targets for 
our June 2004 run. 

2. Within this CMD, we count the number of stars for 
which we expect to obtain photometry down to 1% 
or better, which, according to Fig. 5, will include 
most stars with 14.5 < I < 17. Note that we pre- 
select our targets based on their distance, so that 
stars within this range of apparent magnitude will 
be of spectral type G or later. 

3. As a final step, we perform cuts in V — I color 
to eliminate the redder sequence of background 
evolved disk stars if it is present in the CMD. Stel- 
lar radii of evolved stars are significantly larger than 
their main-sequence counterparts, and thus detect- 
ing planets around evolved stars is virtually impos- 
sible due to the reduced photometric signal depth 
of a transiting planet. We show how we eliminate 
the evolved sequence from consideration in the left 
panel of Fig. 6. 

4. The result of this count approximately corresponds 
to the number of small main-sequence stars we can 
monitor at the 1% photometry level and serves as 
the figure of merit in the cluster selection decision- 
making process. Since the box in the CMD of NGC 
6253 contains more stars (3400) than the one for 
NGC 6134 (2850), NGC 6253 was chosen as our 
observing target for June 2004. The last column 
of Table 1 shows the estimates of the numbers of 
1%-rms stars for our potential target clusters which 
have test data available. 

7. explore/oc observing strategy 

The EXPLORE/OC observing strategy is designed to 
maximize P V i Sl minimize false positives, and to constrain 
field contamination — the issues described in §3 and 4. We 
review aspects of observing strategy that are most impor- 
tant for our project. Some of these are covered in M03 but 
are included here for completeness. We focus in particular 
on considerations necessitated by observing OCs instead 
of Galactic fields. 



7.1. Choice of Filter 

Our photometric monitoring is done in the /-band. The 
shape of a transit in the photometric light curve is de- 
pendent on the filter due to the color dependence of limb 
darkening whose effects are smaller in I than in the bluer 
bands, (see §2 and figure 2 in M03). The transit depth 
is near constant in / when the planet is fully superim- 
posed on the parent star. Because of this "flat-bottomed" 
light curve in /, the shape of the transit makes it easier to 
distinguish planet transits from the signal caused by graz- 
ing binaries (basically a 'pointy' or 'round' eclipse instead 
of a flat-bottomed one) than at bluer bands where limb 
darkening is stronger. Fig. 7 shows a light curve with 
a flat-bottomed eclipse, illustrating that flat bottoms do 
indeed occur at I-band. 

Additional advantages of observing in the /-band are (a) 
increased sensitivity to redder, intrinsically smaller stars 
which offer greater chances of detecting orbiting CEGPs, 
and (b) suffering less extinction due to dust than in the 
bluer bands. 

Disadvantages may include (a) lower CCD quantum ef- 
ficiency in the /-band compared to, e.g., the Z?-band, and 
(b) the occurrence of fringing due to multiple reflections 
and subsequent interference internal to the CCD substrate 
or between the supporting substrate and the silicon. Fring- 
ing is usually more visible in / than in BVR, due to the 
abundant night sky emission lines in the / wavelength 
range. We note that we do not encounter any fringing 
at all with our setup at the Swope Telescope at LCO. 

We also do not change filters during OC monitoring 
since such a strategy would effectively reduce our observ- 
ing cadence (§7.2). 

7.2. Single- Cluster/High- Cadence Observing 

In order to maximize the chance of detecting any ex- 
isting planetary transits, we do not alternate OC targets 
(even though we would increase the number of monitored 
stars that way) but instead observe the same cluster for 
as many hours as possible during the night (Figures 3 and 
4). The main reason for this strategy is to conduct high- 
cadence observing. 

The main goal of this approach is to distinguish a true 
transit light curve from false positives such as grazing 
eclipsing binary stars, an M-star eclipsing a larger star, or 
stellar blends (Seager & Mallen-Ornelas 2003; Charbon- 
neau et al. 2004). Because the total duration of a short- 
period planet transit is typically a few hours, with ingress 
and egress as little as 20 minutes, high-cadence observ- 
ing is essential for well-resolved light curves for a limited- 
duration observing run where only two or three transits 
are expected. A well-resolved light curve with good pho- 
tometric precision can be used to derive astrophysical pa- 
rameters of the planet-star system from the light curve 
alone (e.g., Seager & Mallen-Ornelas 2003) which is useful 
in both ruling out subtle false positives such as blended 
eclipsing binaries, and in obtaining an estimate of planet 
radius. In particular, the density of the parent star is of 
interest for distinguishing between a planetary transit in 
front of a main-sequence star, and the case of a late-type 
dwarf orbiting a giant star. The star's density, however, 
(1) can only be calculated from photometry data alone 
when assuming a stellar mass-radius relation, and (2) is 
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sensitively dependent on the full duration of the transit 
(including ingress and egress), and the duration of totality 
only. 

The flatness of a light curve during the out-of-eclipsc 
stages of a system offers another means of separating plan- 
etary transits from stellar eclipses, as illustrated in Sirko & 
Paczyhski (2003) and Drake (2003). Short-period binary 
stars will have gravitationally distorted, non-spherical 
shapes which will result in a constant sinusoidal brightness 
variation of the light curve with a maximum at quadrature. 

We also do not change targets during the course of an 
observing run of ~20 nights or less (see Fig. 2). The 
justification for this strategy is simple: to maximize P V i S . 
From Panel b of Fig. 2, one can see that the typical val- 
ues for < P vis > (P v i s averaged over all periods between 
1 and 5 days) of a ~ 20-night observing run with some 
holes due to weather will reduce the estimated number of 
detected planets to 50 - 70% of the "theoretical" value as 
calculated in §3.1. Panel c shows that the efficiency, i.e., 
how much is added to < P V i S > per night, will peak at 
around f 8 nights for perfect conditions, justifying our goal 
of observing every cluster for around 20 nights in a row. 

Alternating cluster targets was suggested by Janes 
(1996). Street et al. (2003) adopted an alternating cluster 
strategy, and while their detection algorithm could find 
transits, they found that having only 4 to 6 data points 
observed during transit was a limiting factor in both the 
detection S/N and in discriminating against false positives. 
Furthermore, while alternating cluster targets may provide 
more monitored stars, this strategy will favor only the 1- 
to 2-day period planets if the observing run is not long 
enough (Figure 3). 

Targeting only one cluster during the night further al- 
lows us to keep the stars in the targets OCs on our images 
at exactly the same place on the chip (to within less than 
1 arcsec). This helps us simplify the photometry pipeline. 
In addition, cosmetic problems with the CCD, such as bad 
columns or bad pixels, will eliminate the same stars in ev- 
ery exposure. 

7.3. Dynamic Observing and Optimization of Available 
Telescope Time 

We use a real-time approach to maximizing P vis if the 
allocated observing time is significantly larger than 20 
nights, e.g., > 30 nights, based on detecting a single, full 
transit. 

Panel b of Fig. 2 illustrates that the probability of de- 
tecting an existing single transit (dashed line) will reach 
about 65-70% after around 10 nights of continuous ob- 
serving with 10.8 hours per night. As our data reduction 
pipeline allows us to do practically real-time data reduc- 
tion, we can inspect our highest-quality light curves for the 
existence of a single transit after around 10 nights. If, at 
that point, we do not see any indication of a single transit 
anywhere in our data, we will move on to the next tar- 
get and observe it for the remainder of the allocated time. 
This approach is essentially a comparison of probabilities: 
the probability of detecting two transits in a new cluster in 
the remaining observing time versus the probability (given 
no transits observed so far) of detecting two transits in the 
current cluster if we monitor it for the rest of the available 
observing time. 



7.4. Different Observing Strategies 

In this Section, we describe how different arrangements 
of observing nights affect P V i S . 

At private observatories (such as LCO), different longer- 
term projects requiring many nights may compete for time 
at smaller telescopes such that their allocation of nights 
needs to be split. We explain below how different ways 
of dividing observing time between our project and oth- 
ers affects our likelihood of detecting existing planetary 
transits. 

In Fig. 4, we illustrate the efficiency of a number of 
different observing strategies which may result from such 
split-time arrangements. The solid line in all four panels 
corresponds to P V i S (2 transits detected) of an observing 
run of 20 uninterrupted nights with 10.8 hours of observing 
each night. 

In Panel a, the dotted line corresponds to P vis of an ob- 
serving run spread over 40 nights (10.8 hours per night), 
during which we observe only for the first two nights out 
of every four. P vis is approximately the same as the one 
for 20 consecutive nights. The P V i S averaged over all pe- 
riods (1 day - 5 days), < P V i S >, of the 20-consecutive- 
nights observing run is 0.681. The same < P vis > for the 
2-nights-on, 2-nights-off strategy over 40 nights is 0.666. 
We note that the 2-on, 2-off strategy may impose difficul- 
ties in (1) the period determination due to aliasing effects 
(see below), and (2) the loss of observing time per night 
due to the drift of the sidereal time over the course of such 
a long observing run. 

In Panel b, the dotted line showcases the result of ob- 
serving only the first half of every night for 40 nights in 
a row. The likelihood of detecting existing transits is re- 
duced significantly (< P vis >^ 0.437). For a strategy of 
observing a third of every night for 60 nights, as shown by 
the dotted line in Panel c, < P vis > goes down to 0.007. 

Note that none of these numbers takes into account the 
drift of the sidereal time which would reduce the number 
of hours of observability during the night as a function of 
declination of the target. As a result of the sidereal drift, 

< P V i S > would be reduced from 0.748 for a run of 20 
consecutive nights to 0.705 for a run of 40 nights with the 
2-nights-on, 2-nights-off strategy for NGC 6208, assuming 
it is perfectly centered in RA at the midpoint of the hy- 
pothetical observing run. We calculated similar decreases 
(on the order of 5% or less) in < P vis > when comparing 
the two observing strategies for the other clusters in Table 
1. 

Finally, Panel d illustrates the aliasing effect of only 
observing 2 out of 4 nights. The dotted line in Panel d 
corresponds to the probability of detecting two existing 
transits from which the period can be correctly determined 
when applying the 2-nights-on, 2-nights-off strategy over 
the course of 40 nights. < P vis > of the dotted line is 
0.356, meaning that only about half (0.356/0.666) of all 
transit observations would result in a correct calculation 
of the period, whereas the rest would suffer from aliasing 
effects. Note that this ratio is sensitively dependent on the 
period itself, as illustrated by the dotted line. For compar- 
ison, a strategy of f-night-on, f-night-off would produce 

< P V i S > (two transits observed) of 0.677, but a < P V i S > 
(no aliasing) of only 0.286, meaning that a larger fraction 
(l-(0.286/0.677)~ 58%) of observed transits would result 
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in an incorrect calculation of the period. The strategy of 
continuously observing for 20 nights will give a < P„, s > 
(no aliasing) of 0.406, i.e., an correct estimate of the period 
in 0.406/0.682-60% of the cases. 

We thus conclude that while having 20 consecutive, un- 
interrupted nights is clearly the most favorable solution, 
we can tolerate the strategy where we observe 2 out of ev- 
ery four nights without a significant loss in < P V i S >, but 
which will increase the probability of aliasing effects in the 
period determination. 

7.5. Contamination by Galactic Field Stars 

Estimates of background or foreground stellar contami- 
nation to OCs are valuable since they are the basis upon 
which statistical estimates of planet frequency among OC 
members are based, regardless of whether a planet was de- 
tected or not. In order to get a handle on contamination, 
we observe two control fields per target cluster at the same 
Galactic latitude, approximately a degree away from the 
OC. These observations are ideally taken in V and /, using 
the same exposure time as for the cluster field, and taken 
in the same weather and seeing conditions. To first order, 
the excess number of stars in the cluster field will be rep- 
resentative of the number of cluster members, subject, of 
course, to uncertainty due to fluctuations of background 
and foreground star counts. 

Figures 8 and 9 show this approach for estimating con- 
tamination for the observed OC NGC 2660. Fig. 8 com- 
pares the stellar density (measured in units of stars per 
100 pix x 100 pix on the CCD with 13.0 < I < 17.0) 
as a function of radial distance from the CCD center of 
the cluster image of NGC 2660 (solid line) and two con- 
trol fields (dotted and dashed lines) at the same Galactic 
latitude offset by 1 degree in the sky in either direction. 
The comparison between the CMDs of the cluster and con- 
trol fields is shown in Fig. 9. Although the cluster main 
sequence is not clearly visible in its CMD, one may nev- 
ertheless see a higher density of stars with respect to the 
control field CMDs at colors red-ward of V — I — 1.2, as 
well as a red clump at around I — 13 and V — I — 1.4. The 
total number of stars within 13.0 < I < 17.0 in the cluster 
field is around 3500 stars versus 2700 and 2900 stars in the 
two control fields. The contamination over the entire CCD 
field is thus around 80%, and approximately 30% towards 
the center of the field out to a distance of around 4 arcmin. 

Figures 10 and 11 illustrate how much more severe this 
contamination can be, using NGC 6208 as an example (for 
which we only have data for a single control field). Fig. 
10 compares the stellar density (same units as Fig. 8) 
as a function of radial distance from the CCD center of 
the cluster image of NGC 6208 and a control field at the 
same Galactic latitude offset by 1 degree in the sky. Here, 
the cluster excess stars do not seem to be very centrally 
concentrated (cf. Fig. 8). Finally, the comparison be- 
tween the CMDs of the cluster and control field show a 
slight excess of stars in the cluster CMD at bright mag- 
nitudes (Fig. 11). These excess stars (located around 
/ — 13.0, V — I — 0.7) are evenly distributed over the 
cluster field, and are approaching the bright limit of our 
photometry (see Fig. 5). The total number of stars in the 

13 IRAF is distributed by the National Optical Astronomy Observatories 
in Astronomy, Inc., under cooperative agreement with the NSF. 



above magnitude range in the cluster field is around 6200 
stars versus 6000 stars in the control field. This would 
amount to a contamination of 97% over the entire field, 
and of around 85% in the inner 5 arcmin. This heavy con- 
tamination and the associated high density of the region in 
which NGC 6208 is located was noticed by Lindoff (1972) 
and reiterated in Paunzen & Maitzen (2001). This cluster 
contamination of 97% is similar to the 94% contamination 
of NGC 6819 estimated by Street et al. (2003). Taking 
into account this high rate of contamination, only a few 
hundred stars of high relative photometric precision are 
actually cluster members. Detecting transits around field 
stars, however, is still useful in a number of ways, outlined 
in §1. 

8. explore/oc photometric data reduction 
methods and spectroscopy follow-up 

The EXPLORE/OC strategies concerning photometric 
data reduction and spectroscopy follow-up work are de- 
scribed here in a brief, preliminary way. More detailed de- 
scriptions will follow along with the presentations of our 
results of the individual OCs. 

8.1. Photometry Data Reduction Pipeline 

After the standard IRAF 13 image-processing routines, 
our stellar photometry for the reduction of individual im- 
ages is performed by an algorithm which will be described 
in detail in an upcoming publication (Yee et al. 2004, in 
preparation), and is outlined in principle in §4.3 of M03. 
We will only provide a very brief overview here. 

At the heart of our aperture photometry algorithm is 
the accurate placement of the aperture relative to the cen- 
troid of the star under investigation. This is an important 
issue due to the relative brightness of the sky with respect 
to the monitored stars. To minimize the contribution of 
sky noise and other systematics, we use a relatively small 
aperture (2-3 seeing disks), which further improves pho- 
tometry in the situation of moderate crowding (with star 
separations of a few seeing disks; see §3.3). To achieve 
the accurate placement of the aperture crucial for obtain- 
ing high-precision relative photometry, we use an iterative 
sinc-shifting technique to re-sample every star individu- 
ally such that the central 3x3 pixels are symmetrically 
located about the centroid of the respective star's PSF. 
Performing this shift for every object in the frame is then 
equivalent to using an identical placement of the aperture 
masks for every object, ensuring proper relative photome- 
try. With such re-sampling, aperture photometry of differ- 
ent aperture radii can be performed simply by using inte- 
ger pixel masks of various sizes. Sinc-function re-sampling 
is an ideal method for shifting an image that is Nyquist 
sampled since it preserves resolution, noise characteristics, 
and flux (Hemming 1977; Yee 1988). 

Our relative photometry is then performed by iteratively 
determining the most stable stars within subregions of the 
CCD field. All other stars within the same subregion are 
shifted to the photometric system of these reference stars, 
thereby using iterations to minimize the scatter and to re- 
move outliers from the calculation of the photometric shift. 
The number of iterations, criteria for outlier removal, size 

, which are operated by the Association of Universities for Research 
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of the subregions, and minimum number of stars per sub- 
region are parameters that vary for each dataset. Some 
of the light curves produced by this algorithm are shown 
as examples in Fig. 7 and illustrate our potential to de- 
tect 1% amplitude signals within the intrinsic scatter of 
the high-precision photometry for the target magnitude 
range. 

8.2. Spectral Type Determination Follow-Up 

We determine spectral types for our planet candidate 
stars to provide an independent measure of their sizes 
which may help break degeneracies in the photometric so- 
lution such as period aliasing or stellar blends, and may 
thus determine whether or not costly RV follow-up work 
is desired (Seager & Mallen-Ornelas 2003; Torres et al. 
2004a, b). We have obtained spectral data using a variety 
of instruments which include the Boiler & Chivens Spec- 
trograph and the IMACS Multi-Object Imaging Spectro- 
graph (both on the Magellan 6.5m Telescopes), as well 
as the Wide-Field Re-imaging CCD in Grism/Multi-slit 
mode on the LCO du Pont 2.5m Telescope. We are cur- 
rently analyzing spectral data for our potential candidates 
(examples in Fig. 7) from our work on NGC 2660 and 
NGC 6208 to determine the exact nature of each of the 
systems. Preliminary results are given in Fig. 7. 

Spectral type determination of non-planet-candidate 
stars in the field will give estimates of the foreground red- 
dening along the line of sight and differential reddening 
across the field, and provide an independent check on the 
determination of cluster distance by isochrone fitting. Fur- 
thermore, the knowledge of the spectral types of a repre- 
sentative set of stars (tens or hundreds of stars) will pro- 
vide an additional means of estimating contamination of 
the sample by Galactic field stars and will allow us to de- 
termine the parent sample of non-cluster members. 

9. SUMMARY 

Open clusters are regarded as suitable planet transit 
monitoring targets because they represent a relatively 
large number of coeval stars of the same metallicity located 
at the same distance (§2). Four groups are now monitor- 
ing over a dozen open clusters for short-period transiting 
planets (see §1). 

We reviewed the main challenges facing transit searches 
(§3 and OC surveys in particular §4). In addition to the 
difficulties involved in any transit search, they include: 

• The relatively low number of stars at high rela- 
tive photometric precision (1-1.5%) compared to 
Galactic field surveys of roughly the same magni- 
tude range: ~ 5,000 compared to <~ 50,000 stars re- 
spectively (though the difference in field size is not 
taken into account here). This number is similar to 
the number of stars obtained by the 6x6 degree 2 
shallow transit surveys of brighter stars. 

• The severe contamination by Galactic field stars, 
up to 97% in our clusters for stars at 13 < I < 17. 

• Differential reddening may be problematic in fitting 
isochrones. 



Just like for field transit surveys, OC transit surveys 
need to maximize the number of stars with high photomet- 
ric precision, maximize the probability to detect an exist- 
ing transit, and not be swamped by false positive transit 
signals. 

We presented aspects of the EXPLORE/OC planet 
transit survey design that were considered to meet some 
of the major challenges facing transit surveys. Target se- 
lection is a key aspect to survey design with the number 
of factors over which to optimize (richness, observability, 
age, distance, foreground reddening) actually limiting the 
number of available targets for a given observing time and 
Galactic location. We choose high-cadence observing in or- 
der to sample transits well enough to easily rule out false 
positives such as grazing eclipsing binaries, and to use the 
unique solution method (Seager & Mallen-Ornelas 2003) 
to estimate planet and star parameters. We have shown 
that with an adopted exposure time and a given telescope, 
the distance of the cluster can be chosen to target certain 
spectral types. For the EXPLORE/OC project we do not 
alternate OCs in a given observing run but instead remain 
on one cluster in order to maximize finding planet transits. 
We have shown that this strategy optimizes the probabil- 
ity to detect an existing transiting planet with periods of 
2-5 days days if observing runs are around 20 days. The 
single cluster approach, together with near real-time data 
reduction, allows us to use our dynamic observing strat- 
egy for long observing runs (> 30 nights): if a single full 
transit is not seen within 10 days, the strategy is to move 
on to another cluster. 

EXPLORE/OC is the only OC planet transit survey 
operating in the southern hemisphere. We have presented 
some preliminary data on the OCs NGC 2660 and NGC 
6208 in order to illustrate the main challenges facing clus- 
ter surveys as well as to illustrate our survey design strat- 
egy. Our /-band, high-cadence photometric monitoring 
with the LCO lm Telescope typically attains 1% precision 
in our relative photometry for around 3000-5000 stars per 
OC field in the range 14.5 < 7 < 17 with 5 min exposures. 
For a cluster at a distance of 1 kpc and Eb-v of 0.2, this 
magnitude range corresponds to a range of spectral types 
between mid-G to early M. We have obtained data on three 
additional open clusters: IC 2714, NGC 5316, and NGC 
6253, and plan to target 4-5 more clusters. 

With the ~12 OCs currently being monitored and an- 
alyzed by the four existing OC surveys, there is a good 
chance that some short-period planets will be detected in 
the near future. Because of the potentially large contami- 
nation, and poor availability of physical data on many clus- 
ters in the literature, any detected planets should individu- 
ally be confirmed as cluster members. Furthermore, char- 
acterizing of the cluster parameters is important (Burke 
et al. 2004). With a limited number of stars per clus- 
ter, severe contamination from field stars, and considering 
the finite magnitude range for which high-precision pho- 
tometry can be obtained, only several hundred to a few 
thousand cluster members are monitored with high enough 
photometric precision to detect planet transits; neverthe- 
less, planet transits detected in the contaminating field 
stars are also useful. If one optimizes the important se- 
lection criteria, partly due to the paucity of old clusters, 
most of the suitable OCs for photometric planet searches 
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and radial-velocity follow up can be searched with a rea- 
sonable amount of time and effort. 
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Potential Open Cluster Targets 
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31 


46 


357.9 


-1.6 


+0.07 


7.347 


2 





a Limiting absolute / magnitude to which we can observe with a photometric precision of 1% or better for 300s exposure time at the Swopc 
lm Telescope obtained during photometric conditions and good seeing. This value is obtained by conservatively (cf. §3.1) assuming that the 
apparent Iu m it = 17, and that Aj = 1.94E B -v (Schlegel et al. 1998). 

b Richness class as given in Janes & Adler (1982); Cox (2000) if available. Range: 1 (sparse) to 5 (most populous). Should be regarded as a 
lower limit to the actual richness of the cluster since it depends on the depth of the study from which it was derived (see §7). 

c The approximate number of main- sequence stars (if available) for which we expect to achieve a relative photometric precision of 1% or 
better for 5-min exposures with the Swope Telescope (see §6.3 and Fig. 6 for details). Should be regarded as a lower limit to the number of 
stars around which wc arc able to detect planetary transits (cf. §3.1). 

d Prcviously observed cluster; see von Braun et al. (2004) and Lee et al. (2004) for preliminary results on NGC 2660 and NGC 6208, respectively. 
Wc note that NGC 2660 (our first target) was chosen for its estimated richness and its observability given the allocated observing time. It 
turned out to be a non-optimal target due to its larger distance and correspondingly brighter limiting absolute / magnitude. 



Note. — This table shows our previously observed OCs plus a number of potential target clusters which we chose based on the 
criteria outlined in §5 and §6. Data were taken from the WEBDA database; metallicities from Twarog et al. (1997), available at 
http://obswww.unigc.ch/wcbda/fch_twarog.html. The column headers arc cluster name, distance in parsecs, foreground reddening, limiting 
absolute / magnitude, a, 5, Galactic longitude, Galactic latitude, mctallicity, logarithm of the age (in years), the value for the estimated 
richness class, and the approximate number of stars in the field with relative photometric precision of 1% or better. 
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100 
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Fig. 1. — Depth of transit signal for transiting planets with different radii as a function of MK spectral type and corresponding stellar sizes 
(from Cox 2000) based on geometric arguments only. The diagonal lines indicate the amplitude of the transit signal in the light curve of a 
given planet-star combination. For instance, a Jupiter-sized planet would cause a 0.01 mag dip in the light curve of a GO star, but only a 
0.003 mag dip in the light curve of an AO star. 




period (days) 



Fig. 2. — Probability P v i s of detecting existing transiting planets with different orbital periods. P v i 3 is calculated with the requirement 
that two full transits must be observed. Panel a: P„, s of detecting 2 transits of an existing transiting planet with a period between 1 and 
5 days after 21 (top curve), 14 (second curve from the top) and 7 (bottom curve) consecutive, uninterrupted nights of observing (10.8 hours 
per night). The difficulty of detecting some phases is shown by the dips in the curves (e.g., orbital periods of an integer number of days may 
always feature their transits during the day and are therefore statistically harder to detect). All phases are averaged over for each period. The 
second curve from the bottom shows the real P v i s for our monitoring study of NGC 2660 (19 nights of 7-8 hours per night, with interruptions 
due to weather and telescope scheduling; see Fig. 7). Panel b: The mean P v i a (averaged over 1 day < P < 5 days) as a function of number 
of consecutive nights in an observing run. The solid line is for the requirement to detect two transits and the dashed line for one transit. This 
figure indicates how much the likelihood of finding existing transits grows with an increasing number of nights of observing. Panel c: Run 
efficiency (defined as < P v i s > divided by the number of observing nights) as a function of run length. For the two-transit requirement (solid 
line) and, an observing run of 18 nights is most efficient. For the single transit requirement, the efficiency decreases monotonically with the 
number of nights since additional nights have progressively lower probabilities of detecting "new" transits. 
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Fig. 3. — Pvi-Sphaged as a function of period (in days) for detecting planetary transits in phased data. PiraSp/i aS e<2 is calculated to be 
likelihood that at least N = 20 (solid line), 40 (dotted line), or 60 (dashed line) in-transit points are are accumulated for observing runs of 
different lengths and different observing cadences. The number of points per transit required for a detection is dependent on both SNR and 
exposure time. Panels a and b compare Pvis p h a sed f° r a 21-night (10. 8h) observing run with a cadence of 5 minutes (panel a) and 15 minutes 
(panel b). Panels c (5-minute cadence) and d (15-minute cadence) illustrate the same for a observing run of 40 nights. See text (§3.2) for 
discussion. 
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Fig. 4. — Dependence of P v i a upon period (days) when using different observing strategies. The solid line in all four panels corresponds 
to P v i s (two full transits) of a 20-night (10. 8h) uninterrupted observing run. In Panel a, the dotted line corresponds to P v i 3 of a 40-nights 
observing run during which we only observe the first two out of every four nights. < P v i a > (periods between 1 and 5 days) of the 20-night 
observing run is 0.681. < P v i s > for the 2-nights-on, 2-nights-off strategy over 40 nights is 0.666. The dotted line in Panel b shows < P v i a > 
when observing the first half of every night for 40 nights in a row: < P v i 3 >~ 0.437. When observing a third of every night for 60 nights, as 
shown by the dotted line in Panel c, < P v i s > goes down to 0.007. Finally, Panel d illustrates the aliasing effect of only observing 2 out 
of 4 nights. The dotted line represents the probability of two observed transits being consecutive, as a function of period. Averaged over all 
periods, only about half of all detected pairs of transits would be consecutive when following the 2-nights-on, 2-nights-off strategy. Note that 
these numbers are slight overestimates (few percent) because they do not account for the drift of sidereal time that would affect a specific 
target's observability. For details, see text (§7). 
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Fig. 5. — Photometric precision of night 15 of our monitoring run of NGC 6208. In this diagram, slightly more than 5000 stars have 
photometry of precision 1% or better. This rms is measured as the scatter around the mean magnitude of the star under investigation. The 
1%-photometry stars cover a magnitude range of slightly more than 2.5 mags. The "Z"-shaped feature for stars brighter than / ~ 14 is due 
to the onset of saturation for some of the stars in some of the images in the time series. The clustering of stars around log rms ~ is caused 
by crowding effects when for some of the images, faint stars are blended together with nearby bright stars. 



NGC 6253 NGC 6134 




Fig. 6. — This figure illustrates our approach concerning final target selection, based on our own test data. To choose between the two 
potential target clusters of different distances and foreground reddening estimates, we counted the number of stars in the boxes ("usable 
stars") on the CMDs; the magnitude cuts are conservatively representative of the range within which our monitored stars have photometric 
rms of less than 1% (see Fig. 5). For NGC 6253, we applied a color cut such that our estimate does not include the evolved background 
sequence visible in the CMD to the red side of the box at V — I ~ 1.4, since transiting planets around evolved stars are not detectable due 
to the large radius of the parent stars. For NGC 6134, we do not see an evolved background sequence and thus increased the color range to 
V — I ~ 2.5, in part because the foreground reddening estimate is higher for this cluster (see Table 1). We note that saturation of our stars 
sets in at I ~ 14.5. For details, see §6.3. Estimates for "usable" stars for our other OC targets are given in Table 1. 
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Fig. 7. — Examples of light curves derived from our high-precision relative photometry (see §8) of our NGC 2660 data. Every panel of one 
of the light curves represents the data taken during a single night, starting with night 1 on the bottom left. Night 2's data are shown in the 
panel directly above it, night 3 above that and so on. No data were obtained during nights 13-15 due to telescope scheduling, and nights 6 
and 12 were only partially useful due to weather. All three displayed light curves show the low-amplitude, transit-like signal we are looking 
for in our survey. They are, however, most likely caused by non-planetary phenomena such as a larger-sized companion (left panel) or grazing 
binaries (middle and right panels). Our preliminary work on spectral type determination indicates that star 10099 (left panel) is an early G 
star, star 9079 (middle) is a late A star, and star 13909 (right) is somewhere between F2 and F5. 
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Fig. 8. — This figure compares the stellar density (measured in stars per 100 pix X 100 pix on the CCD with 13.0 < I < 17.0) as a function 
of radial distance from the CCD center of the NGC 2660 open cluster image (solid line) and two control fields (dotted and dashed lines) at 
the same Galactic latitude offset by 1 degree in the sky in either direction in Galactic longitude. The contamination is around 80% over the 
whole field of the CCD, and approximately 30% for the inner ~ 4 arcmin. For details, see §7.5. 




Fig. 9. — This figure shows the CMDs of the field centered on NGC 2660 (left panel) and two control field (middle and right panel) at the 
same Galactic latitude offset by 1 degree in the sky. 
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Fig. 10. — This figure compares the stellar density (measured in stars per 100 X 100 pix on the CCD with 13.0 < I < 17.0) as a function 
of radial distance from the CCD center of the NGC 6208 open cluster image (solid line) and a control field (dotted line) at the same Galactic 
latitude offset by 1 degree in the sky. The contamination is around 97% when integrated over the entire CCD, and around 85% towards the 
inner 5 arcmin. For details, see §7.5. 
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Fig. 11. — This figure shows the CMDs of the field centered on NGC 6208 (left panel) and a control field (right panel) at the same 
Galactic latitude offset by 1 degree in the sky. The comparison between the CMDs shows a slight excess of stars in the cluster CMD at bright 
magnitudes. These excess stars (located around / ~ 13.0, V — I ~ 0.7) are evenly distributed over the cluster field, and arc approaching the 
bright limit of our photometry (see Fig. 5). Previous studies (Lindoff 1972; Paunzen & Maitzen 2001) already mentioned the difficulty in 
separating the cluster main sequence from the Galactic disk population. For details, see §7.5. 



